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The discovery of the top quark in 1995 1 completed 
the third generation of quarks. After years of its dis- 
covery, the top quark plays an important role in theo- 
retical extensions of the standard model (SM). Its large 
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mass gives the top quark a special position within the 
standard model and may shed light on the dynamics of 
electroweak symmetry breaking. One SM extension, top- 
color assisted technicolor [5] , predicts new strong dynam- 
ics which accounts for the spontaneous breaking of elec- 
troweak symmetry and explains the large mass of the top 
quark. This model predicts a vector particle (Z'), which 
couples primarily to the third generation of quarks and 
has no significant couplings to leptons. The existence 
for a narrow {Y z 1 = 0.012 Mz<) width Z' resonance de- 
caying to ti pairs, using leptophobic topcolor model [5] 
as the reference, has been probed since the beginning of 
Tevatron operations both at CDF [J] and DO Other 
theories [6] of physics beyond the SM predict heavy reso- 
nances that add a resonant part to the SM ti production 
mechanism. 

This Letter presents a search for narrow resonant 
states decaying to top-antitop pairs. In the leptopho- 
bic topcolor model, top quarks decay as in the SM via 
the weak interaction, nearly always to a W boson and a 
b quark. W bosons decay into light fermion-antifermion 
pairs: a leptonic decay (32.4%) into a charged lepton and 
a neutrino; or hadronic decay (67.6%) into an up-type 
quark and a down-type quark. All previously reported 
searches have been analyses of top-antitop decays in the 
lepton plus jets channel, where one of the W bosons de- 
cays leptonically (to an electron or a muon) and the other 
W decays hadronically. This channel features a clean sig- 
nature due to the presence of a lepton in the final state, 
and has a branching ratio of 29%. The result presented in 
this Letter is an analysis of ti decays in the all-hadronic 
channel, where both W's decay hadronically. Because 
this topology features only multiple hadronic jets in the 
final state, it is subject to a considerable multijct QCD 
background. We demonstrate that this background can 
be controlled and significantly suppressed with a careful 
event selection. Analysis of ti decays in the all-hadronic 
channel is advantageous for several reasons: the channel 
offers the largest branching ratio (46%) of any of the ti 
final states; there is no unobservable neutrino in the final 
state, which permits improved resolution of the ti invari- 
ant mass; finally, this sample is orthogonal to that of 
previous analyses - the result presented in this Letter is 
complimentary to the previous results in the lepton plus 
jets channel. 

pp collision events analyzed in this paper were pro- 
duced at the Tevatron collider at a center of mass energy 
of 1.96 TeV and were recorded by the CDF II detector [7]. 
The data sample corresponds to a total integrated lumi- 
nosity of 2.8 fb _1 . CDF II is a general purpose particle 
detector. It consists of high precision tracking systems 
for vertex and charged particle track reconstruction, sur- 
rounded by electromagnetic and hadronic calorimeters 
for energy measurement, and muon subsystems outside 
the calorimeter for muon detection. CDF II uses a cylin- 
drical coordinate system with azimuthal angle <fi, polar 
angle 9 measured with respect to the positive z direc- 
tion along the proton beam, and the distance r from 



the beamline. The pseudorapidity, transverse energy and 
momentum are defined as 77 = — In Ptan(f )] , E t — E sin 6 
and P t = Psin#, where E and P are the energy and mo- 
mentum of an incident particle. 

The data were collected using a multijct on line event 
selection (trigger), which is implemented in three stages. 
For triggering purposes, the calorimeter is subdivided 
into a 24 x 24 grid of towers in r/-4> space. At level 
1, we require at least one trigger tower with transverse 
energy E\ ow > 10 GeV. At level 2, we require the sum of 
the transverse energies of all the trigger towers, ^ El° w , 
to be > 175 GeV and the presence of at least four clus- 
ters of trigger towers with Ef s > 15 GeV. Finally, at 
level 3 we require four or more reconstructed jets with 
raw(uncorrected) E t > 10 GeV, where jets are identi- 
fied as clusters of energy depositions in the calorimeter 
using a fixed cone ( AR = ^/(Acj)) 2 + (A77) 2 ) = 0.4) al- 
gorithm [5] . The efficiency of this trigger selection on all- 
hadronic ti events is about 80%. The main background 
present in this data sample is QCD multijet production. 

The jet energies are corrected for calorimeter response, 
multiple interactions and energy radiated outside the jet 
cone [?J]. Jets originating from a b quark are identified 
by the SECVTX [TU] algorithm, which searches for tracks 
with non zero impact parameter that result from the dis- 
placed decay of B hadrons inside the jet, and fits the 
tracks to a common vertex. If this vertex is significantly 
displaced from the primary interaction point, the jet is 
tagged as a 6 jet. 

Events compatible with the signal are selected by re- 
quiring six or seven jets with \rj\ < 2 and corrected 
E t > 15 GeV. To remove leptonic ti decays, we veto 
events with well identified leptons [llj or with signifi- 
cant imbalance in transverse momentum [T2]. After all 
the preselections defined above, the SM ti contribution 
to the data sample is expected to be very small (0.3%). 
To enrich the signal presence in the data sample we have 
to apply additional cuts, which we describe later in the 
paper. 

The distinctive feature of this analysis is the use of like- 
lihoods calculated by integrating signal matrix elements 
both to perform ti invariant mass reconstruction and to 
suppress the overwhelming background. The full expres- 
sion of the likelihood for a given event with jet momenta 
configuration j to be the result of SM ti production and 
decay is given by: 

p( -t, s _ V"^ f dz a dz bf( Z a)f( Z b) TT I" d 3 p, 

nJ\m)-2^J 4Ea E b \v a -v b \ 11 (27r)32^ 

coram 1=1 J 

x \M(m lP )\ 2 (27r) i S^(E F - Ej)TF(j\p)P t (p) 
<Jtot{m)e(m)N combi 

, where z aj b, v a>b and E a ^ are the fractional momenta, ve- 
locities and energies carried by partons a, b and f(z a ,b) 
are the parton distribution functions of colliding proton 
and antiproton, j(p) are jets (partons) four-momenta, m 
is the top quark mass, M(m,p) is the SM ti leading 
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order matrix element, Utot(m)^{'m) is the SM tt produc- 
tion cross section times the selection efficiency both as a 
function of m, and Ep(Ej) is a generic notation for the 
four- vector of the final (initial) state and Pt(p) for the 
transverse momentum of the ti system. The sum is per- 
formed over all jet to parton assignments N com u- The 
probability that a parton with energy E p is observed as 
a jet with energy Ej is given by the transfer function, 
TF(j\p), which is parameterized as a function of parton 
energy and pseudorapidity. Transfer functions are de- 
fined individually for the jets associated with b and light 
jets, as they have different response in the calorimeter. 
We construct the transfer functions from the events that 
have all the jets uniquely matched to each individual par- 
ton within a cone of AR = 0.4. For each energy and i] 
region we use smoothed histogrammed distributions of 
1 — Ejet/Eparton a s TF{j\p) parameterization. A sample 
of fully simulated tt Monte Carlo (MC) events, generated 
using pythia v6.2 [13] with parton showering followed 
by the full simulation of the CDF II detector, and as- 
suming to = 175 GeV/c 2 is used to obtain the Pt(p) and 
TF(J\p) parametrizations. 

The probability density, P(j\m), can be expressed with 
respect to any variable that is a function of parton four- 
momenta, in this case the invariant mass of tt pairs, Ma- 
Integrating Eq. ([T]) over the parton variables times a delta 
function 5{x — Ma{p)) we obtain the probability density 
function for each event once the jets are measured. We 
use the mean value of this distribution as an estimator 
for M a . 

To discriminate between SM background and Z' signal 
events, we calculate event quantities which are sensitive 
to the presence of a signal and use them as inputs for 
a neural network which is trained to separate the signal 
and the background. Keeping in mind that SM tt is one 
of the background samples for Z' events, here we will 
refer to SM tt as the signal sample for the event selec- 
tion purpose only. We train the neural network to select 
events with the presence of tt pairs and to veto dominant 
QCD multijet production. Using SM tt events as signal 
events to optimize the event selection and enrich the tt 
content of the sample accomplishes reasonable results for 
Z' events as shown later in the paper. In addition, this 
choice makes the search unbiased to a specific mass and 
the model of Z' hypothesis used. 

A first set of 10 kinematic variables, summarized in 
Table |TJ has already been shown to be effective [14] in 
reducing the QCD background. Significant distinguish- 
ing features of tt production in comparison to the QCD 
background are high Et jets, dijet resonances from W 
decay and trijet resonances from t decay. The central- 

ity is C = ^j-r. > where v s is the invariant mass of the 

multijet system. The aplanarity is defined as A = §Qi, 
where Q\ is the smallest of the three normalized eigenval- 
ues of the sphericity tensor, M ab = J2j PjPji calculated 
in the center of mass system of all jets, where a and b re- 
fer to the spatial components of the jet four-momentum 



Pj. In Table ni 9* is a jet emission direction, represented 
by the angle between the jet direction, measured in the 
center of mass frame of all jets, and the proton beam 
axis. For the last variable, MED (matrix element dis- 
criminant), we exploit the broad set of information from 
the event about its production and decay through the 
SM it matrix element. For each event we calculate 'the 
minus log probability' of Eq. ([T]) at 9 different top mass 
points, Tfif — 155, 160. ..195 GeV/c 2 , and use their sum 
as the final discriminator. 



TABLE I. Neural network input variables. 



Variable 


Description 




Scalar sum of all jet Et 




As above, excluding two highest Et jets 


c 


Centrality, defined in text 


A 


Aplanarity, defined in text 


M 2 7 n 


Minimum dijet invariant mass 


M max 


Maximum dijet invariant mass 


M 3 7 n 


Minimum trijet invariant mass 




Maximum trijet invariant mass 




Et sin 2 9* for the highest Et jet 


(E* T ) 


Geometric mean of Et of remaining N — 2 jets 


MED 


Constructed from matrix element 



Having defined the variables, to separate tt from back- 
ground events, we use them as inputs to a neural net- 
work [15] with two hidden layers and one output node. 
The neural network is trained on samples of signal and 
background events with 6 < Aj c t s < 7. To model the sig- 
nal events we use the pythia MC generator at leading 
order(LO) to produce SM tt events assuming a top quark 
mass of TOtop = 175 GeV/c 2 and the theoretical cross sec- 
tion of 6.7 pb [15] . We use the multijet data events as 
the background sample for training the neutral net since 
the tt contribution is expected to be negligible. After the 
training the value of the output node, NN out , is used as 
a discriminator between signal and background. Its dis- 
tribution is shown in Fig. [I] In addition, we show the 
comparison of the QCD dominated data, MC generated 
SM tt and MC generated Z' events for one of the input 
variables ^2 Et ■ 

In the final event selection we require a cut on the 
neural net output, NN out >0.93, and at least one jet 
tagged as having originated from a b quark. The neu- 
ral net requirement was optimized to suppress the QCD 
background while enhancing the content of tt events by 
maximizing the SM tt significance. Table [TT] shows the 
selection efficiencies for SM tt and Z' events after final 
event selection cut. There are 2086 events surviving these 
final selection criteria including 680 SM tt events as es- 
timated from the simulated event sample and assuming 
the NLO theoretical cross section [IB]. The remaining 
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Neurei net output 



FIG. 1. ^2 and neural net output value for QCD domi- 
nated data (black points), SM tt (dark grey histogram), QCD 
prediction (light grey histogram) and 700 GeV/c 2 tt resonance 
normalized to the expectation (dashed histogram). 



events are from QCD multijet processes plus a potential 
signal contribution from Z' events. 



the final selected sample. To test our background model, 
we consider several control regions, denned by the neural 
net output value NN out < 0.25, 0.25 < NN out < 0.75, 
0.75 < iViV out < 0.93. For all the regions we find a 
very good agreement between the model and the observed 
data. Figure [2] shows the distributions of ^2 Et and M t j 
in the control region 0.75 < NN out < 0.93. 



i 



$NDF=25.2/45 
p-ob=0.94 



100 200 300 400 500 600 700 80 
lEj [GeV] 




x)NDF = 25.4/39 
prob=0.85 



300 400 500 600 700 800 900 1000 
M t f [GeV/c 2 ] 



FIG. 2. Distributions of ^ Et and M t i in the control region 
0.75 < iViVout < 0.93. 



TABLE II. Table of cross sections, a, and acceptances, e ± 
<5e(tot.), for Z' and SM tt events. 



M z /(GeV/c 2 ) 


a[pb] 


e 




SM tt 


6.7 


3.8 


± 


0.5 


450 


8.96 


4.2 


± 


0.5 


500 


5.66 


4.7 


± 


0.5 


550 


3.40 


5.3 


± 


0.5 


600 


2.09 


5.7 


± 


0.5 


650 


1.31 


5.8 


± 


0.4 


700 


0.78 


5.6 


± 


0.4 


750 


0.47 


5.2 


± 


0.3 


800 


0.28 


4.6 


± 


0.3 


850 


0.16 


4.0 


± 


0.2 


900 


0.10 


3.6 


± 


0.2 



The dominant background is multijet production via 
QCD, where one of the b jets can originate from heavy 
flavor (b or c) quark pair production or from misidenti- 
fied light flavor quark jets. Due to the large theoretical 
uncertainties on the production cross section, we use a 
data-driven approach to estimate the QCD background. 
From a data sample with 4 or 5 jets, which is overwhelm- 
ingly from QCD production (SM tt fraction less than 
5- 10~ 4 ), we build a tag rate matrix. In this procedure, we 
parametrize the probability for each jet to be identified 
as a b jet. The parametrization includes the dependence 
on the transverse energy of the jet, the number of tracks 
associated to the jet, and the number of reconstructed 
collision vertices in the event. Once we define the prob- 
ability for a single jet to be tagged, we can use the tag 
rate matrix to estimate the probability for an event to 
have one or two 6-tagged jets. The tag rate matrix is 
applied to 6 or 7 jet data events before the 6-tagging re- 
quirement to predict the QCD background for events in 



To measure and set the confidence level intervals on 
resonant tt production given the observed M t i spectrum, 
we start with the following likelihood, 

..rii 

L{n\cr,V) = J{e-^ (2) 

i 

which is the prior probability of observing ft, where rii is 
the observed number of events in the ith bin of the M t i 
distribution, and fa is the expected number of events in 
the same bin and is given by /x, = er s (^4 s — A c ° nt )T % s + 
a u (A tt - A% nt )Tl t + N QCD T l QCD , which depends on the 
assumed signal cross section, a s , the SM tt cross section, 
an, the signal and SM tt effective acceptances, A S ,A U , 
the number of expected QCD events, Nqcd, and the 
fraction of events in ith M t i bin, T l s , T£ t , Tq CD for the 
signal, SM tt and QCD distributions, respectively. As 
the QCD background prediction is performed using the 
data sample itself, the presence of SM tt and assumed 
signal events must be subtracted from QCD background 
estimation. We calculate the number of residual con- 
tamination events in QCD background prediction, using 
assumed and theoretical cross sections for signal and SM 
tt events, and their effective acceptances for the residual 
contamination terms A c s ont and A1° nt by applying the tag 
rate matrix to the signal and SM tt samples before the 
tagging requirement, and 

We use Bayes' theorem to connect the likelihood of the 
measurement to the posterior probability density, which 
is used to set the upper limits. 

p{a\n) = J dvp(a,v\n) = J dVL{n\a, v)/p{n){i) 

where 7r(<7, u) is the prior probability density, and p(ii) = 
J dv J daL(a, v\n)Tr(a, V). 
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There are two types of uncertainties we have to con- 
sider. The first type does not change the shape of the M t j 
distribution but only the acceptances. The second type 
affects both the shape and normalization; we'll refer as 
shape uncertainties. Uncertainties that do not change the 
shape of the templates (distributions) are incorporated as 
nuisance parameters and integrated over in Eq. @. In 
this respect, Eq. (|3| includes not only the statistical un- 
certainty of the data, but also the source of systematic 
uncertainties on: signal and SM ti acceptances, SM ti 
cross section, QCD normalization and integrated lumi- 
nosity. 

One of the shape uncertainties we consider is the jet 
energy scale corrections. After the jet energy corrections 
we are left with an uncertainty on the jet energy scale. 
A change in the jet energy scale modifies both the accep- 
tances and the template shapes. To account for shape un- 
certainties and jet energy scale in particular, we generate 
a set of pseudoexperiments using the shifted templates 
and acceptances. This results in a shifted reconstructed 
cross section with respect to the nominal one. The map- 
ping of this shift versus the input cross section provides 
an evaluation of the impact of the jet energy scale un- 
certainty at any given cross section. The complete list 
of the shape uncertainties sources we considered are jet 
energy scale, initial and the final state radiation, and un- 
certainty on proton and antiproton parton distribution 
functions. Assuming that the nature of shape uncertain- 
ties follow normal distribution, we convolute the poste- 
rior probabilities with a Gaussian whose width is equal 
to the quadrature sum of individual shape uncertainties. 

After including the shape uncertainties, the posterior 
density function is used to define the upper and lower 
limits at any given confidence level (C.L.). If the lower 
limit is zero then the data is considered consistent with 
the SM at that level of confidence. We also extract the 
reconstructed cross section as the most probable value of 
the posterior. To obtain the sensitivity of the reconstruc- 
tion algorithm we generate 1000 simulated experiments 
in the signal null hypothesis and extract the 95% C.L. 
expected upper limit, defined as the median of the up- 
per limits distribution. This entire exercise is repeated 
for various resonance masses from 450 GeV/c 2 to 900 
GeV/c 2 . Together with the theoretical cross section ver- 
sus mass curves these limits are used to exclude certain 
mass ranges. 

In this analysis we consider the data gathered by CDF 
II between 2002-2008 at the Tevatron. The M ti distri- 
bution for the 2086 events surviving the final event se- 
lection criteria is shown in Fig. [3] and is consistent with 
the SM expectations. The resulting 95% C.L. upper lim- 
its on a(pp — > Z') ■ BR{Z' — > ti) as a function of M t j 
are shown in Fig. [4] together with expected limits de- 
rived from pseudoexperiments that include the SM back- 
ground hypothesis only. These limits can be used to ex- 
clude a leptophobic topcolor resonance candidates with 
a mass less than 805 GeV/c 2 at 95% C.L., assuming the 
width of the resonance is Tz< — 0.012 Mz>- The previous 



searches were performed in the lepton plus jets channel 
only, and the most recent results were conducted by CDF 
II [T7] and DO [TB]. Using Tevatron data corresponding 
to 1 fb^ 1 and 0.9 fb _1 integrated luminosity respectively, 
they found no evidence for ti resonant production. For 
the same benchmark model of leptophobic topcolor Z' , 
the upper limits were set at 720 GeV/c 2 and 700 GeV/c 2 
for CDF II and DO, respectively. 




MJGeV/c] 

FIG. 3. Reconstructed M t t vs the SM expectation in the 
search region above the 400 GeV/c 2 cut. 






450 500 550 600 650 700 750 800 850 900 

M z - [GeV/c 2 ] 



FIG. 4. Expected and observed upper limits on leptophobic 
topcolor Z' in 2.8 fb -1 of CDF II data. The blue line is the 
median expected upper limit with the assumption of no signal, 
the red line is the observed limit and the black line is the cross 
section prediction for leptophobic topcolor Z' production. 

In conclusion, we have performed a search for a heavy 
resonance decaying into ti using data with 2.8 fb _1 inte- 
grated luminosity in the all-jets channel. No evidence is 



observed and we set upper limits on the production cross 
section times branching ratio at the 95% C.L. For one 
leptophobic topcolor production mechanism, we exclude 
masses up to 805 GeV/c 2 . 
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